Pion and Kaon Condensation at Finite Temperature and Density 
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In this paper, we study 0(2iV)-symmetric </> 4 -theory at finite temperature and density using the 
2PI-1/7V expansion. As specific examples, we consider pion condensation at finite isospin chemical 
potential and kaon condensation at finite chemical potential for hyper charge and isospin charge. We 
calculate the phase diagrams and the quasiparticle masses for pions and kaons in the large- N limit. 
It is shown that the effective potential and the gap equation can be renormalized by using local 
counterterms for the coupling constant and mass parameter, which are independent of temperature 
and chemical potentials. 
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I. INTRODUCTION 

0(iV)-symmetric field theories are playing a signifi- 
cant role as effective quantum field theories in condensed 
matter and high-energy physics. For example, the non- 
relativistic 0(2)-symmetric model is used to describe 
cold atomic Bose gases in the context of Bose-Einstein 
condensation in harmonic traps [l[ and the homogeneous 
dilute Bose gase [2j. Similarly, the 0(4) (non)linear 
sigma model is used as a low-energy effective theory for 
QCD. In massless two-flavor QCD, the global symmetry 
group is SU(2)l x SU(2)r and one then takes advantage 
of the isomorphism between this group and 0(4). This 
effective theory describes the three pions and the sigma 
particle. Similarly, the SU (2) l X Si7(2)fl-syinmetric the- 
ory is a toy model for the description of kaon conden- 
sates [3, 0] in the color-flavor locked phase [f| of high- 
density QCD. 

The thermodynamic functions for hot field theories 
can be calculated as a series in the the coupling con- 
stant by resumming certain classes of diagrams from all 
orders of perturbation theory. However, it turns out that 
the expansion of these functions in powers of the cou- 
pling constant breaks down at temperature unless the 
coupling is tiny [|| 0, [|[ . In order to obtain more stable 
results, one has developed and applied more advanced 
resummation techniques, such as screened perturbation 
theory Q, hard-thermal-loop resummation [To| and re- 
summations based on the two-particle irreducible (2PI) 
effective action [III EH EH EH- These approaches all 
improve the converge of the perturbative series of the 
thermodynamic quantities. 

At large chemical potential and zero temperature, re- 
summation of classes of diagrams is also required [H, Efl 
to obtain the thermodynamic functions as a series ex- 
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pansion in the coupling. However, the perturbative se- 
ries may be better behaved at zero temperature than 
at high temperature as the plasmon effects are differ- 
ent [17] . In this case, and also at nonzero but low 
temperature (T < fi), hard-dense- loop resummation has 
proven to be an effective way of reorganizing the calcu- 
lations fig . Very recently, Ipp et al pj| have presented 
a new method for calculating the pressure in deconfmed 
QCD at all values of the temperature and the quark 
chemical potential, which reduces to the known pertur- 
bative results in the appropriate limits. 

The 1/N expansion in conjunction with e.g. the 1PI 
or 2PI effective-action formalism is another nonpertur- 
bative resummation scheme which has been applied at 
finite temperature [H HH HH HI • A leading-order cal- 
culation typically gives rise to a gap equation for a mass 
parameter, which is straightforward to solve. At next- 
to-leading order (NLO), the resulting equations become 
harder to solve. Due to nontrivial renormalization is- 
sues, the NLO problem using the 1PI l/A^-expansion, 
has been solved rigorously only very recently [22] . In 
particular, it was shown that the thermodynamic poten- 
tial can be renormalized in a temperature-independent 
way only at the minimum. In the case of the 2PI for- 
malism there are additional complications since the gap 
equation for the propagator becomes nonlocal and non- 
trivial to solve. However, the 2PI formalism is superior 
when it comes to the description of nonequilibruim phe- 
nomena [2lj . 

Renormalization of expressions obtained in resumma- 
tion schemes such as the 2PI effective-action formal- 
ism is a nontrivial issue. In these approximations, one 
generally resums only certain classes of diagrams from 
all orders of perturbation theory. Since only a subset 
of diagrams from a given order in peturbation theory 
is taken into account, standard proofs of perturbative 
renormalizability do not apply and one is not guaranteed 
that ultraviolet divergences can be systematically elimi- 
nated. Moreover, it is not obvious that these divergences 



are independent of temperature. Signficant progress on 
these issues has been made in recent years starting with 
the papers by van Hees and Knoll [2J|. They showed 
that the equation for the two-point function in scalar 
<^ 4 -theory can be renormalized by introducing a finite 
number of counterterms. A more systematic study of 
the renormalization issues in scalar </> 4 -theory was car- 
ried out by Blaizot, Iancu and Reinosa [25j]. In partic- 
ular, they focused on the elimination of temperature- 
dependent divergences. This is sufficient to calculate 
physical quantities such as the thermodynamic pressure, 
but only recently has it been shown by Berges et al [26[ 
how to fix all the counterterms needed to calculate the 
proper vertices which are encoded in the effective action. 
These counterterms are local and they are independent 
of temperature and chemical potential. 

The idea of a Bose-Einstein condensate of pions or 
kaons in the core of compact stars has a long history 
and the interest in such condensates is due to their far- 
reaching consequences for these objects. For example, 
the presence of a condensate results in enhanced neu- 
trino cooling (27[. Moreover, in contrast to hadronic 
matter in heavy-ion collisions, bulk matter in compact 
stars must (on average) be electrically neutral and so a 
neutrality constraint must be imposed. 

In the present paper, we discuss pion and kaon con- 
densation at finite temperature and density using the 
2PI-1/7V expansion. Pion condensation and the phase 
diagram of two-flavor QCD have been investigated us- 
ing chiralperturbation theory [2H [29l . [30l [3lj . lattice 
QCD (Hi!, ladder QCD [||, chiral quark model [H[, 
and Nambu Jona-Lasinio models in the mean-field ap- 
proximation [36], [I?], [Hj]. In Ref. [38f , the effects of 
imposing charge neutrality were studied at zero tem- 
perature. There have also been some applications of 
the linear sigma model at finite temperature and den- 
sity in either the large-iV approximation [37j or the 
Hartree approximation [39j. In both cases, the quasi- 
particle masses were calculated, but renormalization is- 
sues were not discussed. In the present paper, we 
generalize part of the finite-temperature calculations of 
Lenaghan and Rischke [23| by including chemical po- 
tentials. More specifically, we study pion and kaon con- 
densation at finite temperature and density using the 
2PI 1 /iV-expansion paying attention to renormalization 
issues. 

The paper is organized as follows. In Sec. II, we 
briefly discuss interacting Bose gases at finite density. 
In Sec. Ill, we consider pions at finite temperature and 
isospin chemical potential. We calculate the phase dia- 
gram and the quasiparticle masses. In Sec. IV, we in- 
vestigate kaon condensation in the presence of chemical 
potentials for isospin and hypercharge. In Sec. V, we 
summarize and conclude. 



II. INTERACTING BOSE GAS 

The Euclidean Lagrangian for a Bose gas with N 
species of massive charged scalars is 



c : uyi>; no,,*, > - ^=(*i + * f i) + ™ 2 *!*i 



'2N 



(1) 



where i = 1,2, ...,7V and $i = (<p2i-i + i<p2i)/V2 is 
a complex field. If H = 0, the theory described by 
Eq. |T]) has (27V — 1)N conserved charges which equals 
the number of generators of the group 0(2N). If H is 
nonzero, the 0(2N) symmetry is explicitly broken down 
to 0(2N — 1). A gas of TV species of bosons can be 
characterized by the expectation values of the different 
conserved charges in addition to the temperature. For 
each conserved charge Qi, one may introduce a nonzero 
chemical potential fii. However, it is possible to specify 
the expectation values of different charges only if they 
commute. The maximum number of commuting charges 
is N [401 and these can be chosen as 



where the the current densities are 



n 



(2) 



(3) 



The incorporation of a conserved charge Qt is done by 
making the substitution 



<9o$; ->■ (do - Hi) $i 



(4) 
(5) 



in the Lagrangian (JTJ) . From the path-integral represen- 
tation of the free energy density T 



0VF = V<5>*V<5> ie -fo dT J d3xC 



(6) 



the expression for the charge density can be written as 



Qt 



(7) 



If we introduce k chemical potentials, the full symmetry 
group is broken down to [0(2)] fc x 0(2N-2k). Urn 2 < 0, 
the 0(2N) symmetry is spontaneously broken down to 
0(2N — 1). Even if m 2 > 0, the symmetry may be 
broken if one of the chemical potentials, fii, is larger than 
a critical chemical potential /i c = 771. In that case, the 
(O|02i— 1|0) 7^ and the corresponding 0(2) symmetry 
is broken. 
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III. PIONS 

In this section, we discuss pions at finite chemical 
potential pi for (the third component of the) isospin 
charge, which corresponds to the current density j M ~ 
</>2<9 M 3 - fod^fo, and N = 2 in Eq. Q. In the vac- 
uum, the 0(4) symmetry is broken down to 0(3) by a 
nonzero vacuum expectation value </>o for <p\. The sym- 
metry is broken either spontanously in the chiral limit 
or explicitly if H ^ 0. 



We first introduce a nonzero expectation value po for 
<f>2 to allow for a charged pion condensate. The field <f>i 
is then written as 



$i = — ={(j} + ip + <f>i , (8) 



where <pi and cf>2 are quantum fluctuating fields. The 
inverse tree-level propagator Dq (u> n ,p) reads 



D Q 1 (uj n ,p) 



( Lu^ + p 2 + m\ 

$<t>aPo 






T!<t>oPa 



+ V ■+ 
2p I uj r , 







-Zpi^n 






(9) 



Ldl+p 2 + ml 



where p = |p|, u n = 2nnT is the nth Matsubara fre- 
quency, and the treel-level masses are 



jr[0 o ,p o ,£>] 

Spo 



= 



(16) 



3A 2 *^ 2 

2iV 00 + 2N Po 



"Mi 



771, = — 



Mi 



A 

2iV~' 
A 

2~/v' 



3A 
2~/V 
_A_ 
2~/V 



2iV 00 + 2N Po 



(10) 

(11) 
(12) 
(13) 



Note that the mass parameter m 2 is negative in the re- 
mainder of this section, as we want to break chiral sym- 
metry spontaneously in the chiral limit in the absence 
of a chemical potential. 

The 2PI effective potential can be written as 



T[<p ,p ,D] 



1 2 

-m 



Pi) 



Po) 



-TjM/Po - Hc t>0 



877 (l 

- -TrlnD- 1 
2 



-TvD^D + $[£>] , 



(14) 



where D is the exact propagator and &[D] is the sum of 
all two-particle irreducible vacuum diagrams. The trace 
is over space-time as well as field indices. The expecta- 
tion values 0o and po satisfy the stationarity conditions 



$r[<f>a,p ,D] 



= 0, 



(15) 



The exact propagator D satisfies the variational equa- 
tion 



ST[(/) ,p ,D} 
SD 



. 



(17) 



Using Eq. (jTTJ) and II = D 1 — D 1 , where II is the exact 
self-energy, the variational equation can be written as 



11(D) = 2 



S$[D} 
8D 



(18) 



The diagrams that contribute to the self-energy matrix 
H(ll!,p) are obtained by cutting a line in the vacuum 
diagrams. 

In Ref. [2lj], the diagrams contributing to &[D] were 
classified according to which order in the 1/N expan- 
sion they contribute. The various terms can be writ- 
ten in terms of O(N) invariants such as Tr(D") and 
Tr(</>gZ) n ), where n is an integer. At leading order in 
1/N the diagrams that contribute are shown in Fig. [TJ 
This diagram can be written as an 0(./V)-invariant term, 
$ LO ~ (TrD) . The vertex gives a factor of 1/N, while 
each trace yields a factor of N. It is therefore a leading- 
order contribution. 

The diagram that contributes to the self-energy II to 
leading in 1/N is shown in Fig. O The particle prop- 
agating in the loop is one of the 27V — 3 particles with 
tree- level masses given by Eq. (JT5J). 
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FIG. 1: Feynman diagram contributing at leading order in 
l/N to 




FIG. 2: Feynman diagram contributing at leading order in 
1/iV to the self-energy Tl(uj,p). 



Their inverse propagator to leading order in 1 /N can 
be written as 



D-\uj n ,p) = P 2 



(19) 



where P 2 



p 2 and nLo(-D) is the tadpole in 



Fig. ^ calculated with the full propagator D. Since 
IIlo(-O) is independent of the external momentum, the 
propagator can be parametrized as 



D 1 (uj n ,p) 



P 2 + M 



(20) 



where M is a constant. The leading order self energy 
can now be written as 



n 



LO 



= A 



1 



M 2 



(21) 



Using the expression for IIlo and the two expres- 
sions (fT9|) and l|20|) for the inverse propagator, we obtain 
a gap equation for M: 



M 



2 = m 2 A + \ 



1 



Q 2 + M 2 ■ 



The sum-integral is defined by 
VA 2 



JQ 



-in 



q =2-KnT 



d d q 



(22) 



(23) 



where d = 3 — 2e and A is the renormalization scale 
associated with dimensional regularization. The integral 
over 3-momentum will be calculated using dimensional 
regularization. Inserting the parametrization (|20[) of D 



into the 2PI effective potential (|14[) and differentiating 
with respect to 4>o and p$, the equations (|15p and (|16p 
can be written as 



= mVo-H + ^o(^o+P6) 



-X<j> 



M 2 



(24) 



= (m 2 -/i?)p + — Po(0o+Po) 



-Ap 



(25) 



f Q Q 2 + M 2 ■ 

Combining the gap equation ([22| and Eq. ([24| gives 

i? = M 2 o . (26) 

Similarly, combining the gap equation (|22| and Eq. (|25|) 
gives 



(M 2 - M 2 )po 







(27) 



We will study the various phases in the following sub- 
sections. 

We next discuss the renormalization of the gap equa- 
tion ([22]) . which can be done by an iterative proce- 
dure given in Ref. [25| (see also Ref. |ll|). By writing 
M 2 = m 2 + IIlo and expanding the right-hand side of 
Eq. (f2"2")l in powers of IIlo, it is clear that it rcsums the 
perturbative daisy and superdaisy diagrams: 



M 



= A 



1 



q Q 2 + m\ 
W 



1 



+n LO A 



4 

JQ 



n LO 

IIlo A 
1 



(Q 2 



ml) 2 



(Q 2 + m 2 ) 



2 A3 



(28) 




FIG. 3: Diagrammatic interpretation of the gap equa- 
tion ([221) for n LO- 



This is illustrated diagrammatically in Fig [3] Renor- 
malization is carried by making the substitutions 



m 2 + >J &m 2 n , 

71=1 

oo 

A + SX n > 



(29) 
(30) 
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where 5m 2 and SX n are counterterms of order A" and 
A n+1 , respectively. Similarly, we write the leading order 
self-energy in a power series in A: 



n LO = n 



LO,l 



n 



LO,2 



(31) 



where the subscript indicates the power of A. Insert- 
ing Eqs. (|29|) - ([3Tj) into equation (j28|) , we can determine 
nLO,n, fir^n and SX n by iteration. The first iteration of 
the gap equation yields 



IIlo,i 



Q 2 



X 2 , SX 1 



Pi) 



(32) 



where we have used that IIlo 



M 2 



The next 



step is to calculate the sum- integral in Eq. (j3"2")l . The 
sum over Matsubara frequencies can be expressed as a 
contour integral in the complex energy plane. By per- 
forming this integral, we obtain 



1 



d d q 



1 



+ M 2 J q {2ir) d + m 2 



1 



= VQ2 +M 2 /T _ 



(33) 



The first term is divergent in the ultraviolet. Calculating 
it with dimensional regularization and expanding in e 
through order e , we obtain 



-finite terms . 



The divergences are removed by the counterterms 



(35) 



Sm 1 = 



m 2 X 
167r 2 e 



16tt 2 6 



(36) 



This renormalization procedure can be carried out iter- 
atively to all orders in A and this determines the coun- 
terterms Sm 2 an d 6X n for n > 1. The counterterms <5tti 2 
and SX n are [25| 



1 



Q 2 



- M 2 
1 

'2^2 



A/ 2 
16tt 2 
dq q 2 



A2_ 

M 2 



1 

- + 1 

e 

1 



o vV + M 2 eV V+I^/T 



(34) 



5m„ = 



8X n = 



A n+1 



(37) 



(16^ 2 e)" ' (16^ 2 e) n ' 

Summing the counterterms to all orders, Eqs. (f2"9")l 
and (l30l) can be written as 



The logarithimic divergence in the sum-integral in 
Eq. (|2"2")) shows up as a pole in e, while the quadratic 
divergence is set to zero in dimensional regularization. 
Inserting the sum-integral (|34| into Eq. (|32)) . we obtain 



n 



A? 



LO,l 



167r 2 e 



Sm 2 H 

1 2N 



Pi) 



16tt ;i € 



167r 3 e 



(38) 



Returning to the nonperturbative gap equation (|2"2"|) . it 
can be rewritten as 



M 2 
A 



m 2 1 

m 2 1 
~ + 2N 



pI) 



pI) 



i 



Q W 
Af 2 
16tt 2 



+ Af 2 
AT 2 



l 

- + 1 

e 



2tt 2 



dq q 2 



o V? 2 + M 2 eV^+^/i 



(39) 



Making the substitutions 



the renormalized gap equation becomes 



M 



X , n, AA/ 2 / A 2 



A 

~2tt 2 



dq q 2 



o y/q 2 + M 2 e vV^/r - l 



(40) 



We finally discuss the renormalization of the effective 
potential in the broken phase. In the large- iV limit, $[-0] 
is given by 



M 2 



(41) 



Using the expression (|4ip and inserting the gap equa- 
tion (f22ll into the expression for the effective poten- 
tial (fT4")l . we obtain 



r[0 Q ,Po,£] = ^ (Af 4 - m 4 ) - _ ^ 



-NM- 



M 



(42) 



The first sum-integral is given by Eq. (|34|) . while the 
second can be found by integrating the first with respect 
to M 2 and reads 



ln(P 2 



M ) = - 



M 4 


(m 2 ) 


'1 


3" 


32^ 2 




e 


2 



T 

^2 



dqg 2 ln(l-e-V?+^/ T ) .(43) 



Inserting the expressions for the sum-integrals into ([421 
and making the substitutions (|38p , we obtain the renor 
malized effective potential 



T[<f> ,p ,D] 



N 
2A 



(M 4 



1 



In 



NM 4 

nm 2 r°° 

lo 



2tt 2 

f CO 



7 2 

1 

M 2 + 2 



vV + M 2 gV^+M^/T _ ! 
~J"dq^hi(l-e-^ r ^ S / T ) , (44) 

where we have omitted a temperature-independent di- 
vergence proportional to m 4 . This corresponds to 
adding a vacuum counterterm AS to the effective po- 
tential. We also note that the effective potential can 
be renormalized in a similar manner in the symmetric 
phase. 



A. Chiral limit 

In this subsection, we discuss the solutions to the gap 
equations and the corresponding phases in the chiral 
limit, i.e. when H = 0. 

First assume that the charged pion condensate po is 
nonzero. Equation (|27[) then immediately gives 



M 2 = p 2 . 



(45) 



Equation 
vanishes: 



Using Eqs. 



26)) then implies that the chiral condensate 

0o = 0. (46) 
(14"5"[) and (|4"6"| , the pion condensate reads 

1 



™tf- m *)-2N 



Q 2 + i4 



(47) 



The pole in e on the right-hand side of this equation 
is removed by making the substitutions (|38p in analogy 
with the renormalization of the gap equation ([2"2")) . After 
renormalization, the gap equation (|47[) can be written as 



pI 



2N 2 



M 

8tt 2 



N 

"^2 



dq q 2 



In — 5- 
1 



1 



(48) 



As the temperature T increases, Eq. (|4"51) shows that po 
goes to zero continuously and so the phase transition 
is second order. The critical temperature T c is given by 
the temperature for which the right-hand side of Eq. (|48|) 
vanishes. In the weak-coupling limit, where A <C 1, we 
can obtain an analytic expression for T c . In that case, 
pi can be neglected in the sum-integral in Eq. (|47p and 
its value is NT 2 / 6. This yields the critical temperature 



Tl = f (M? - m 2 ) 



(49) 



which is in agreement with the calculations of Ka- 
pusta 

The dispersion relations are obtained by analytic con- 
tinuation to Minkowski space, Lo n — > iuj, and then solv- 
ing the equation det D(oj,p) = 0. For T < T c , the dis- 
persion relations are 



J 



(50) 



G 



W2,3(P) 



f + 2^2 + 1 ( M | 



Ml) ± J (Mi- Mir + 8/i?(2/i? 



M| 



M|) 



•4/*fp 



2„2 



(51) 
(52) 



where M 2 = to 2 +I1lo- Note that the a and 7r° are iden- 
tified with the original fields (f>i and (f>4, and therefore 
with u)\(p) and W4(p). On the other hand, the charged 
pions tt~ and n + are given by linear combinations of 
4>2 and 03 and are identified by u)2(p) and ^(p). In 



the minimum of the effective potential (fl4|) . we have 
M| = 2(p 2 I -rh 2 ) and Mf = 0, where to 2 = m 2 + IT 
The dispersion relations in Eq. (|5"lj) then reduce to 



1LO- 



^2,3(P) 



3/x? - to 2 ± J 



2^2 



(53) 



From this equation, we see that the ir + is a massless 
mode, which in the long-wavelength limit behaves as 



u<s(p) = 




— P 



(54) 



This is consistent with the fact that the original 0(2) x 
0(2) symmetry is broken down to 0(2) by the pion con- 
densate. The fact that the Goldstone mode is behaving 
linearly with p makes the system superfluid as the Lan- 
dau criterion is satisfied [431 ] . 

In the normal phase, where 4>a = Po = 0, the disper- 
sion relations are 



^2,3(P) 
U} 4 (p) 



VP 2 + M 2 , 
VP 2 + M 2 ± fij 



(55) 
(56) 
(57) 



vV + M 2 , 

where M is determined selfconsistently by solving the 



gap equation (f40|) evaluated at O = Po = 0: 



M 2 = to 2 



2tt 2 



XM 2 A A? 
16tt 2 V "m 2 
dqq 2 



o vV + M 2 eV^W/T 



(58) 



B. Physical point 

We next discuss the case where H ^ 0. The solution 
to Eq. dHJ) yields 



H 

14 



(59) 



Note that (f>o is nonzero and independent of the temper- 
ature. Assuming a nonzero pion condensate, the mass 
M again is given by Eq. (|45]). Using Eqs. (J45j) and ([59]) . 
the gap equation (|25|) can be written as 



2iV 



if 2 

A 4 / 



io 2 + 



-2 -(so) 



The gap equation is renormalized in the same way as 
the gap equation in the chiral limit. 

The expressions for the dispersion relations are very 
involved and they are not particularly illuminating. 
We therefore restrict ourselves to list the quasiparti- 
cle masses. These are found by solving the equation 
detD(uj,p = 0) = and read 



W 2 ,3(P = 0) = M 2 , 3 = 



\ 



7 2 

2^ 



oji(p = 0) = Mi = 



±1^(5^ -2m 2 ) 2 -12^ 



'N v 2 ' 

u>i{p = o) = m 4 = m , 



(61) 

(62) 
(63) 



where m 2 — m 2 + IIlo- The ir + is again a massless dispersion relations around p = 0, it can be shown that 
Goldstone mode. By expanding the equation for the this mode is linear in p and hence the system is super- 
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fluid. This mode reflects the breaking of the residual 
0(2) symmetry and is a conventional Goldstone mode. 

We next discuss the system in the normal phase, 
where po = 0. Eliminating 0o from the gap equa- 
tion J2H), one finds 



2N 



M d -m z ) = 



El 

Af 4 



2N 



1 



M 2 



(64) 



This equation is renormalized in the usual manner. The 
renormalized gap equation determines the mass M as a 
function of temperature. We see that the equations (|6T)|) 
and (|64f match onto each other continously at the phase 
boundary, where M — fii as they should. The dispersion 
relations in the symmetric phase 



wa,a(p) = VP 2 + M 2 ± h , 
w 4 (p) = VP 2 + M 2 , 



(65) 

(66) 
(67) 

Note that all four modes are massive as the 0(2) sym- 
metry has been restored. 

C. Numerical results 

We next discuss the determination of the parameters 
in the Lagrangian (JTJ. In the vacuum, we have po = 
Hi = and f„ = (f>o- At the tree-level, we have 

(68) 
(69) 

H = Uml . (70) 

The mass parameter to 2 and the coupling constant A 
can then be expressed in terms of the observable masses 
of the sigma and the pion: 



"i ■ 

l 2 ' 



to 



A 



N 



(m 2 a - ml) 



2 ) 



J TV 



(71) 
(72) 



Eqs. ([70 |) -(|72 ]1 can now be used to determine the con- 
stants iJ, to 2 , and A, given the values m a = 600MeV, 
m^ = 139MeV,and / w = 93MeV. In the chiral limit, 
we find to 2 = -(300MeV) 2 and A = 83.2466. At the 
physical point, we obtain to 2 = — (230.5MeV) 2 and 
A = 78.7788. Once we take into account quantum fluc- 
tuations, Eqs. and ([55)) are modified. After renor- 
malization, we find 
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16tt 2 
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(73) 
(74) 



Note that the renormalization scale A is a free parameter 
and in the chiral limit, it does not enter our expressions. 
At the physical point we chose A 2 = to 2 /e, so that the 
tree-level relations still hold without having to change 
the parameters m 2 and A. 

In Fig. 31 we show the medium dependent meson 
masses at the physical point as functions of the isospin 
chemical potential at zero temperature. For /i/ > \i c = 
139 MeV, they are given by Eqs. ![oT ]) -([g3 ]) . while for 
lii < fi c , they are given by Eqs. (f65"l) - (f6"7|) with p = 0. 
For fii > [x c , the mass of 7r + vanishes. 
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FIG. 4: Pion and sigma masses at the physical point as 
functions of the isospin chemical potential fii an< l T = 0. 



In Fig. [5j we show the pion condensate as a function 
of temperature normalized to the condensate at T = 0. 
In both the chiral limit and at the physical point, the 
pion condensate breaks an 0(2) symmetry. From uni- 
versality arguments, wc know this transition is of second 
order, which can be seen from Fig. (O. The critical ex- 
ponent v = 1/2, which is expected from a mean-field 
calculation. See Ref. [44| for a NLO calculation of the 
critical exponents using the 2PI formalism. 

In Fig. [SI we show the phase diagram for pion con- 
densation. The solid line corresponds to the chiral limit, 
while the dashed line corresponds to the physical point. 
The two curves meet at large \ii since the temperature- 
independent value of H becomes irrelevant. 

In Fig. we show the meson masses as functions of 
temperature for fixed isospin chemical potential fii — 
200 MeV. For T < T c , they are given by Eqs. ([6T]) - 
(j63j) . while for T > T c , they are given by Eqs. (|S5 > (|6? )) 
with p = 0. The critical temperature is T c — 223 MeV. 
Note that one of the pions is massless below the critical 
temperature 
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FIG. 5: Charged pion condensate as a function of temper- 
ature normalized to the condensate at T = 0. The upper 
curve is the chiral limit and the lower curve corresponds to 
the physical point. 
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FIG. 6: Phase diagram for pion condensation. The upper 
curve is the chiral limit and the lower curve corresponds to 
the physical point. 



IV. KAONS 

In the color-flavor locked phase of high baryon-density 
QCD, the lightest excitations are the four charged and 



neutral kaons, i"T , K° and K°. The low-energy ef- 
fective theory is described by the Lagrangian {T]) with 
H = and N = 2 fs]. In this section, we discuss kaon 
condensation in the presence of chemical potentials fly 
for hypercharge Y and /i/ for isospin. The kaons can 
be written as a complex doublet, (<&i,$2) = {K°,K + ). 
We then have the following relation between the various 
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FIG. 7: Meson masses as functions of temperature at the 
physical point. The isospin chemical potential is fii = 200 
MeV. T c = 223 MeV. 



chemical potentials: 



1 

2 



(M2 - Ml) , 
■ 



(75) 
(76) 



We first write the complex fields <&i and $2 in terms of 
the charged and neutral kaon condensates </>o and po and 
quantum fluctuating fields 4>\ — 4>4- 



$1 = i (0o + 0i + ifo) , (77) 



$2 = — ^ (po + 03 + i<t>i) 



(78) 
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The tree-level propagator is 



D 1 (uj n ,p) 
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(79) 



where the tree-level masses are 
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(80) 
(81) 
(82) 
(83) 
(84) 



Note that in the remainder of this section, the mass 
parameter m 2 is positive. The 2PI effective potential is 
given by 



f 



T[(j)Q,p,D] = ^ (m 2 - p{ 



2^2 



1/2 2\ 2 

- (TO - /x 2 j p 



- — (0o + P 2 ) 2 + -Trln.tr 1 
8A^ ; 2 



(85) 



To leading order in N, the contribution to <&[£>] is again 
given by the diagrams in Fig. [TJ The stationarity condi- 
tions (fl"5 j) - ([T6|) for the effective potential are 



-Xp 



M 2 



(87) 



The self-energy to leading order in 1/N is again diago- 
nal and contributing diagram is shown in Fig. O The 
expression is given by Eq. (|22|) . Using Eq. (|86|) for the 
minimum of the effective potential, we obtain 



M 



Mr 



(88) 



We first assume that p\ = p 2 . In this case, the La- 
grangian (QJ has an extended SU{2) x U(l) symmetry 
instead of 0(2) x 0(2). If p\ > m 2 this symmetry is 
broken down to U(l). We have complete symmetry be- 
tween (j)Q and po and it is then possible to rotate e.g. the 
po condensate away. If pi > p 2 (to be specific), po = 
and 4>o ^ only if pf > to 2 . In this case the 0(2) x 0(2) 
symmetry is spontaneously broken down to 0(2) by the 
condensate. 

In analogy with Eq. (f2"4")l , the kaon condensate (j>o sat- 
isfies the gap equation 



2N 



Mr 



TO 2 ) 
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M 2 



(89) 



(to 2 - p 2 )^ + {4>l + pl) 



+A0O 



2 + M 2 ' 



A 



= (m - p 2 )p + —po 



Pi) 



The phase transition is again of second order. In the 
weak-coupling limit, an analytic expression for T c analo- 
(86) gous to Eq. (|4"9"1) is obtained by the substitution pi — > pi . 

In the broken phase, the dispersion relations can be 
written as 



^1,2 (p) = 

^50) = 



p 2 + 2p\ + ~(M 2 + M 2 ) ± ^{M 2 - M 2 ) 2 + 8p 2 (2p 2 + M 2 + M 2 ) + Ap\p 2 



P 2 + Mi ± M2 , 



p 2 + p\ 



(90) 

(91) 
(92) 
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where M 2 = mf + IIlo The neutral kaons K° and 
K° are identified with the linear combinations of the 
fields 4>\ and </>2, and therefore with ui\{jp) and uj^ip). 
Similarly, the charged kaons K~ and K + are given by 
linear combinations of <fo and <f>4 and are identified with 
u> 3 (p) and u>i(p). 

At the minimum of the effective potential, we have 
M 2 = 2{n\ — in 2 ) and Mf = 0, where in 2 = m 2 + IIlo- 
Eq. ([50]) then reduce to 

wi,2(p) = y^p 2 + 3/i? - m 2 ± y ;3/xf - to 2 ) 2 + 4p 2 // 2 . 

(93) 

This shows that W2(p) and thus if is a massless exci- 
tation which is linear for small p (cf. ([51)1 ) and has the 
quantum numbers of if . In the case where //i ^= U2, 
this is the only Goldstone boson which corresponds to 
the broken 0(2) symmetry. In the case where /ii = [12, 
wa{p) is another massless modes,which is quadratic for 
small p, which has the quantum numbers of K + . Thus 
there are two massless modes despite the fact that there 
are 3 broken generators [H, [4(| [47| ■ The splitting be- 
tween K + and K~ is due to the fact that the chemical 
potentials break the discrete symmetries C, CP, and 
CPT [47|. This result is in agreement with the counting 
rules derived by Nielsen and Chadha [Hj], which state 
that the massless modes with a quadratic dispersion re- 
lation must be counted twice. Finally, we notice that 
only in the case where ^ ^2, the condensate implies 
that the system is superfluid. Whenever the Goldstone 
mode with a quadratic dispersion relation is present, the 
Landau criterion for superfluidity can never be satis- 
fied 0,113. 

We next discuss the system in the normal phase, 
where 4>q = 0. The dispersion relations reduce to 



wi, 2 (p) = VP 2 + M 2 ± mi , (94) 
w 3 ,-i(p) - vV + M 2 ± M2 , (95) 
w B (p) - vV + M 2 , (96) 

where M 2 again satisfies the gap equation (f58|) . Since 
the symmetry is restored, all excitations are gapped. 

In Fig. [3 we show the quasiparticle masses as a func- 
tion of fix at zero temperature. These are given by 
Eqs. ([90 |) -(|92 p with p = in the broken phase and by 
Eqs. (f94 |) - (|96|) with p = in the symmetric phase. In the 
broken phase, the neutral kaon K° is a massless Gold- 
stone mode. We have set fi2 = A*i/2 and chosen the 
renormalization scale A = M in the symmetric phase. 
The phase diagram is the same as in the pion case with 
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FIG. 8: Kaon masses at zero temperature as functions of 
the chemical potential fii, where fii > p,2- 

V. SUMMARY 



In the present paper, we have calculated to leading 
order in the 2PI-1/7V expansion finite density and finite- 
temperature effects in interacting Bose systems with 
pion and kaon condensation. In particular, we have ob- 
tained the phase diagrams and quasiparticle masses. De- 
pending on the magnitude of the chemical potentials, the 
number and nature of the Goldstone modes vary. In all 
cases, our results are in agreement with general theorems 
on the occurence of massless excitations [45j, |48|, |49(. 
This is in contrast with e.g. the Hartree approximation, 
where the Goldstone theorem is not satisfied. Further- 
more, we have shown that is is possible to renormalizc 
the theory in a medium-independent way by eliminating 
the divergences through a redefinition of the mass and 
coupling constant. 

The present work can be extended by calculating the 
1/N corrections either in the context of the standard ef- 
fective action or by using the 2PI effective-action formal- 
ism. Using the conventional 1PI effective action, the cal- 
culations would be similar to those recently carried out 
in Refs. [2l],[22| ■ Calculations within the 2PI effective ac- 
tion will be more difficult as the resulting equations form 
a set of coupled integro-differential equations. These 
equations must be solved numerically, which is nontriv- 
ial. However, progress has been made recently [14| in 
massive scalar field theory with vanishing chemical po- 
tential. 
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